Calorimetry of Bose-Einstein condensates 
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We outline a practical scheme for measuring the thermodynamic properties of a Bose-Einstein condensate as 
a function of internal energy. We propose using Bragg scattering and controlled trap manipulations to impart 
a precise amount of energy to a near zero temperature condensate. After thermalisation the temperature can 
be measured using standard techniques to determine the state equation T(U,N,uj). Our analysis accounts for 
interaction effects and the excitation of constants of motion which restrict the energy available for thermalisation. 



I. INTRODUCTION 

The harmonically trapped Bose-Einstein condensate (BEC) 
system is well-isolated from its environment and the thermal 
state can be characterized by the parameters of total atom 
number (N), internal energy {U) and trap potential frequency 
(uj). To date the internal energy of a Bose gas has been im- 
practical to measure experimentally with a useful degree of 
accuracy. On the other hand, when a discernible thermal frac- 
tion is present the temperature can be quite accurately deter- 
mined by absorption imaging after expansion [1]. For this rea- 
son the temperature dependence of BEC thermodynamics, e.g. 
condensate fraction versus temperature, is quite well known, 
whereas the energy dependence has hardly been studied. The 
few experimental studies conducted have suffered from large 
uncertainties such that any form of quantitative comparison 
with theory was not possible [2]. More recent studies have 
used an increase in temperature to signify the imparting of 
energy 0], but have not considered the quantitative rela- 
tionship between these two quantities. 

Knowledge of the energy dependence of BEC thermody- 
namics is of wide spread interest. The A-transition in He was 
so named according to the peculiar shape in the specific heat 
capacity of the system [5]. Additionally a detailed description 
of energy dependence would be useful for discriminating be- 
tween finite temperature theories of ultra-cold Bose gases. For 
the case of degenerate Fermi gases a heat capacity measure- 
ment has been made by the Duke group [6] by manipulating 
the trap potential in a manner similar to what we consider here 
for Bose gases. Additionally, recent work by the Heidelberg 
group [q\ has examined a precise method for measuring the 
temperature, and used this to confirm the deviation of the heat 
capacity of a Bose gas from that of a classical gas for a con- 
stant background noise source. However, their input heating 
rate was unknown, so that they were unable to quantify the 
heat capacity. The MIT group H |^] have measured heating 
in a BEC by stirring it with a blue-duned focused light beam, 
and used those measurements to distinguish between theories 
for drag forces. 

Here we propose using two mechanisms for transferring 
a precise amount of energy to a BEC at near zero temper- 
ature to establish the relationship between energy and tem- 
perature. For making calorimetric measurements, one ideally 
would like a well-defined reservoir to transfer heat to the sys- 



tem of interest. The isolation of ultra-cold atom experiments 
make such an approach impractical, however an irreversible 
work process, such as that done by a spinning paddle wheel 
in a fluid, is a convenient method for transferring energy into 
these systems without changing the external constraints. 

Our main concern in this proposal is to develop and analyse 
precise ways of imparting energy to the system, and to charac- 
terize the portion of this energy that is irreversible. While our 
analysis here focusses upon Bragg scattering, and expansion 
from a trap, one could envisage doing this with other methods, 
e.g. general perturbations of the trapping potential or stirring 
with a focused light field. The main requirement is that the 
energy transfer is accurately calculable. 



II. PRECISE ENERGY TRANSFER 

A. Overview of proposal 

In the next subsections we analyze two methods for pre- 
cisely transferring energy (ftrans) to the system to be rether- 
malized. In Sec. Ill CI we consider the use of Bragg scattering, 
and in Sec. Ill Dl we consider the sudden expansion from a har- 
monic trap. Having added this energy, the final temperature 
of the system (after it has returned to equilibrium) can be ac- 
curately measured. Knowledge of the irreversible work done 
on the system and the rethermalized temperature establishes 
an equation of state relationship of the form T(U), where T 
is the temperature and U is the internal energy (relative to the 
energy of the T = ground state). An important considera- 
tion in equating U to £trans, is that £trans must only consist of 
the irreversible work done on the system. In particular, energy 
transferred to the Kohn mode (of harmonically trapped gases) 
needs to be excluded (as we discuss below). 



B. Zero temperature formalism 

We consider our initial system to be a Bose gas at zero tem- 
perature, where the condensate is essentially pure. The con- 
densate orbital satisfies the time-independent Gross-Pi taevskii 
equation 
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the harmonic trapping potential, and {uii ,^2,^3} are the trap 
frequencies along the coordinate directions. The quantity \x is 
the chemical potential, N the number of particles in the ini- 
tial (pure) condensate, and Uq = ^ah 2 /m is the interaction 
strength, where a is the s-wave scattering length. 

The energy of the ground state is given by the energy func- 
tional 
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Many-body corrections to the ground state will in general be 
important, however our interest here lies in understanding how 
much energy is transferred to the system rather than the abso- 
lute energy, for which the energy functional will suffice. 

On several occasions we will have cause to make use of the 
Thomas-Fermi approximate solution to the Gross-Pitaevskii 
equation ([B 
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with the chemical potential determined by 
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where ui = (U1W2W3) 1 / 3 (e.g. see [10]). This approximation 
is found by neglecting the kinetic term in the Gross-Pitaevskii 
equation, which is usually much smaller than the potential 
and interaction contributions. Comparisons with experiments 
have shown the Thomas-Fermi approximation to be a good 
description of T = condensates (improving in accuracy as 
the number of particles in the condensate increases, e.g. see 
Ml 

The Thomas-Fermi ground state has an energy of 
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calculated from Eq. (0, with the kinetic term neglected. 



C. Method I: Energy Transfer by Bragg scattering 

Our first scheme for imparting energy uses the well- 
understood process of Bragg scattering 111 lL [l3ll that is 
routinely used in labs for manipulating BECs. We consider 
the situation where a Bragg pulse is used to first-order Bragg 
scatter a fraction a of the condensate (at rest) to momentum 
state Kb (where b is the reciprocal lattice vector of the Bragg 
potential). Note that the harmonic trap remains on during the 
Bragg scattering and the subsequent dynamics of the system 



as it rethermalizes. We also assume that the duration and in- 
tensity of the Bragg potential are chosen so that all other or- 
ders of scattering can be neglected, yet the scattering can be 
considered approximately instantaneous on the timescale of 
condensate evolution (e.g. see Ref. 111411 ). The matter wave 
field at the conclusion of the scattering is then given by 



*i(x) = vT^J 9 (x) + V7ie ibx * s (x), 



(7) 



where we have assumed that the size of the condensate is large 
compared to 1/b, so that the wave-packets centered at mo- 
menta and Kb are orthogonal. 



1. Kohn Mode 

The initial state (Q has a momentum expectation of = 
ahb per particle and according to Kohn's theorem 111511 this 
will lead to an undamped dipole oscillation in the harmonic 
external potential. Transforming to the (non-inertial) time- 
dependent centre-of-mass frame of reference, this oscillation 
can be removed, while leaving the Hamiltonian for the sys- 
tem unchanged (see lfl5lo . In making this transformation an 
amount of energy corresponding to the energy of the dipole 
oscillation, i.e. 



E D = Np 2 /2m = a 2 NKoj b , 



(8) 



is removed, where u>b = Kb 2 /2m is the Bragg recoil fre- 
quency. Because this energy is locked into the centre of mass 
oscillation it is not available for rethermalization. 

However, if the motion of the scattered atoms enters into 
a sufficiently anharmonic region of the trapping potential, 
then the Kohn mode energy will be available to rethermal- 
ize. We expect this to be strongly dependent on the manner in 
which the harmonic potential is made, but this effect should 
be clearly observable as a decay in the centre-of-mass (COM) 
oscillation of the system. 



2. Transferred energy 

The energy transferred to the condensate in the lab frame 
by the Bragg scattering is calculated in Appendix I A II After 
subtracting the energy locked into the Kohn mode we obtain 
that transferred energy available for rethermalization is 
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where £ ; ° t = (N 2 U /2) J d 3 x |* ff (x)| 4 is the interaction en- 
ergy of the ground state. In the term proportional to huib in 
Eq. (O, the a contribution arises from the transfer of kinetic 
energy to the particles, where as the —a 2 part accounts for 
the energy locked into the Kohn mode. The term proportional 
to 2£° lt describes the additional energy arising from interac- 
tions due to the creation of a coherent density fluctuation in 
the system. We have taken this term to lowest order in the 
small parameter X/L, where A = 27r/6 is the wavelength of 
the density fluctuation and L is the size of the condensate. For 
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6.6 x l(r 29 Jf« 4.8^K 


4.3 x 1(T 31 J^ 0.031ajK 


87 Rb 


9.8 x 1(T 30 J^ 0.71/iK 


7.3 x l(T 31 Jf« 0.053^K 



Table I: Energy scales per particle for typical experimental parame- 
ters. Each case considers a 10 6 atom condensate in an isotropic 50 
Hz harmonic trap. We have taken the usual scattering lengths for 
each atom llOfl and have taken b to be that for counter-propagating 
light fields of wavelength A = 589 nm and A = 789nm for the 
sodium and rubidium cases respectively (i.e. b — 4tt/X). 



the trap centre may also change if the trap frequencies of the 
final trap are different to those of the initial trap. For this rea- 
son we do not give explicit expressions for the dipole energy 
here. 



2. Transferred energy 



Using the results of the Appendix I A 21 we find that the en- 
ergy transferred and available for thermalisation is 



more details we refer the reader to Appendix I A II The size 
of these kinetic and interaction contributions to the energy are 
compared in Table [I] We note that for both cases considered 
the interaction contribution is < 7% the kinetic contribution. 
So for many cases ignoring the interaction term in Eq. (0 will 
be a good first approximation. 

We note that the maximum transfer of energy occurs when 
a = 1/2, since the energy available for rethermalization is 
proportional to the momentum spread in the initial state dT) 
which is maximized for 50% scattering. For a > 1/2 the en- 
ergy transferred by the Bragg scattering is increasingly locked 
into the Kohn mode. 



iV/iTF 



2 - 5 7 6/5 



(12) 



where jj = uij/uij, 7 = -^717273 and /itf is the Thomas- 
Fermi chemical potential of the initial condensate. Unlike the 
Bragg case, the transferred energy is unbounded since the con- 
densate can be allowed to expand for arbitrarily long time pe- 
riods. 



3. Analytic solution 



D. Method II: Energy Transfer by Expansion from trap 

Our second method of energy transfer is to suddenly turn 
off the harmonic trapping potential (of initial frequencies 
{ojj}) for a period of time t on , allowing the condensate to ex- 
pand, before the potential is reinstated (with final frequencies 

To quantify the energy transfer in this process we need to 
consider the condensate dynamics after trap release. Castin 
and Dum have shown that such a system will undergo a self- 
similar expansion 11611 . In particular, if the trapping frequen- 
cies are time dependent, then the condensate density remains 
as a Thomas-Fermi profile, but with time-dependent widths 
that evolve according to 



where the equation of motion for the Xj are 
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1. Kohn Mode 

When the trap is turned off at t = the condensate begins to 
expand and it will also fall a distance d = a g t 2 /2 and acquire 
a velocity v = a g t, where a g is the acceleration due to gravity. 
At time i on , when the trap is restored, d and v will manifest 
themselves as energy locked into the dipole mode fl34ll . We 
note that, depending on how the harmonic trap is produced, 



While the energy transfer generally requires us to solve the 
ordinary differential equations ( fTTT ). an approximate solution 
exists for the case of an elongated (cigar) trap. Here the pa- 
rameters Aj^ and X z specify the system at time t. Defining 
e = uji/oj z and t = ui±t we can find an approximate solu- 
tion (vg\ 
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X z = 1 + e 2 [r tan -1 r- In v / TTr 5 ] +0(e 4 ). (14) 

Thus the energy available for thermalisation given the same 
initial and final trap parameters is 
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With T on = 10 I t on . 



III. APPLICATION OF CALORIMETRY 

We now consider the application of our calorimetry scheme 
to an ideal trapped Bose gas with critical temperature given by 
T c = || [A7C(3)] 1/3 . For T < T c , we have that the energy 
of this gas is given by 
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A. Bragg limitations 

In order for the Bragg scheme to be capable of probing up to 
the transition region, we require that maximum energy trans- 
fer is greater than the energy content of the gas at the critical 
point (i.e. £trans(a = 1/2) > U(T C )). Evaluating this in- 
equality sets the following constraint on the system size and 
trap 
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Thus for N or u> too large, first order Bragg will not provide 
sufficient energy to heat the BEC from T = to T = T c . For 
example, 150 x 10 3 rubidium-87 atoms in a 25 Hz trap is about 
at this threshold. For sodium-23, the higher kinetic energy of 
the Bragg scattering (hu>b) allows larger/tighter systems to be 
used. One could circumvent this limitation by using second 
order Bragg scattering to impart a larger amount of energy. 

We note that given the high recoil temperatures (relative 
to usual T c values) associated with the light used to Bragg 
scatter it may seem surprising that it is not always possible 
to Bragg scatter sufficient energy to take the system to the 
critical point. However, it should be kept in mind that the 
recoil temperature is defined for an atom with a mean energy 
equal to the recoil energy in each degree of freedom, whereas 
for our Bragg scattering process this energy must be shared 
between the 6 degrees of freedom of the atom (i.e. 3 kinetic 
and 3 potential), and at most only 50% of the total energy 
transferred is available for rethermalization. 



B. Error analysis of Bragg energy transfer 

Error in the transferred energy will likely be dominated by 
shot-to-shot variation in the number of atoms in the initial con- 
densate and error in the fraction of scattered atoms. Here we 
denote these errors as AN and Aa respectively and account 
for their effect on our ability to know the precise amount of 
energy transferred to the system. Linearizing Eq. (O, we find 
that the energy transferred is affected by these quantities ac- 
cording to 
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where we have neglected the contributions of the interaction 
term since it is typically an order of magnitude smaller. In 
Fig. Q]we consider the effect of these errors in the use of our 
scheme to determine the heat capacity. 

Additional error will arise because the initial condensate 
temperature is non-zero, however since the energy is a rapidly 
increasing function of T [35], this error is typically several 
orders of magnitude smaller than the errors due to uncertainty 
in atom number and scattered fraction. 

In Fig. [T]we also compare U(T) for the interacting ver- 
sus ideal system. The interacting properties are calculated 
using the Hartree-Fock-Bogoliubov theory in the Popov ap- 
proximation (e.g. see HI [ll, Hi HI HHj). We see that for 



=? 3 

CM 

i 

O 

>, 

E? 
<5 

LU 1 



(a) / 






// 






// 






































































20 40 
T [nK] 



Figure 1: (a) Internal energy (U) versus temperature, (b) The spe- 
cific heat capacity (AU/ AT). Ideal gas (solid line), Hartree-Fock 
Bogoliubov-Popov calculation of interacting gas (dashed). Energy 
corresponding to a = 0.5 is shown as horizontal dotted line in (a) 
and grey region indicates the uncertainty in the transferred energy for 
the ideal case given that there is a 5% uncertainty in the total number 
of atoms (i.e. AN = 0.05 x N), the fraction scattered is accurate to 
the 5% level (i.e. Aa = 0.05 x a). System parameters: 1.5 x 10 5 
Rb atoms in an isotropic 25 Hz trap. 



T < T c ~ 60 nK the dependence of U on T for the inter- 
acting system is noticeably distinguishable from the ideal gas 
(i.e. beyond the limits of the error in energy transfer). This 
suggests that interaction effects could be experimentally mea- 
sured using this technique if sufficiently good reproducibility 
of initial condensate number and Bragg scattering precision 
can be obtained. 



C. Error analysis of energy transfer by trap expansion 



As discussed in section UlD 2l imparting energy by trap ex- 
pansion has no upper bound to the amount of transferred en- 
ergy. We would expect that the shot to shot variation in atom 
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numbers will dominate the error budget, however there may 
be other important considerations relating to the particular 
way the trap is produced. 

For the analytic case given in Sec. Ill D 31 we examine the 
sensitivity of the energy transfered to variations in atom num- 
ber (AN) and errors in the axial (Aoj z ) and radial (Auj±) trap- 
ping frequencies. By linearizing Eq. ( TT3T > we find 
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Assuming that the dominant uncertainty is atom number vari- 
ation, we see that the relative uncertainty in transferred en- 
ergy is roughly proportional to the relative uncertainty in atom 
number. For the general case given in Eq. (fT2b . the error asso- 
ciated with trap uncertainties requires a numerical integration 
of the differential equation (fTTT i. However, the (likely domi- 
nant ) error associated with number uncertainty is still given 
by the first term in Eq. ( fT9l ). 



IV. CONCLUSIONS 
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These integrals are all of order X/L where L is the spatial 
size of the condensate and A = 2ir/b is the wavelength of 
the Bragg induced density modulation. In experiments typical 
values are X/L ~ 1/50, so these approximations are well 
satisfied and higher order terms can be ignored. 
With the above approximations we obtain 
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are the single-particle and interaction energies of the ground 
state respectively, with £, — c " 
energy is given by 



£g p + £P t . Thus the transferred 
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The Thomas-Fermi approximation to the ground state (01 
gives 



We have presented two practical schemes for performing 
calorimetry on a Bose-Einstein condensate system. It is clear 
from our results that reasonably accurate calorimetry mea- 
surements could be made using Bragg scattering or by con- 
trolled expansion from a confining potential. We have charac- 
terized the sensitivity of these methods to typical experimen- 
tal uncertainties in atom number and have also shown that it 
should be feasible to measure interaction effects on the ther- 
mal properties of a Bose gas. 

Our scheme also presents a rather well-defined initial con- 
dition for studying non-equilibrium dynamics. There is con- 
siderable interest in the dynamics of the thermalisation, as we 
expect there will be a crossover from coherent to incoherent 
dynamics. This topic is of significant current interest (e.g. see 
SllllIllll|23,|2i,|2i,[M[Il) andwill be the subject of 
future work using classical field methods (e.g. 132113311 ). 



Appendix A: ENERGY CALCULATIONS 

1. Bragg transferred energy 

Here we calculate the total energy transferred to the con- 
densate by Bragg scattering in the lab frame. This is found 
by evaluating Straus — ^[^i] — ■E'l^g]- Assuming that the 
original and scattered wave packets are well separated in mo- 
mentum space, we can make the approximations 
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£P t = -iV MTP . (A8) 



2. Trap expansion transferred energy 

Having solved Eqs. ( fTTT i. the condensate density is then 
given by 



|*(x,t)| 2 = 
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a. Expansion 

We envisage turning the trap off suddenly. In the Thomas- 
Fermi approximation, for the ground state, the kinetic energy 
is negligible meaning that the system only has interaction en- 
ergy given by the generalized form of Eq. dA6l ) 
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The condensate begins to expand: the interaction energy is 
converted to kinetic energy (Skin)- The energy balance is de- 
termined by £kin(i) = £int(0) — £int(X)> an d using the time- 
dependent density of Eq. (|A91 > we find 
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Now consider turning the trap back on suddenly at time 
t = t on , with some possibly different trap frequencies {w'j}- 
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Defining the ratio of new to old trapping frequencies as jj = 
ujj/uij, the potential energy at the instant the new potential is 
turned on is given by 
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and the total energy of the system will be the sum of the in- 
teraction energy at t = and the potential energy at t — t on , 
i.e. 
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The using the Thomas Fermi result in Eq. (O we find the 
approximate ground state energy in the final trapping potential 
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where we have defined 7 as the geometric mean of the 
7y . Thus the energy added by the sudden release and re- 
application of the harmonic trapping potential is 
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